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Abstract 

The electromagnetic calorimeter at the BaBar detector, part of the asymmetric B Factory at SLAC, 
measures photons in the energy range from 20 MeV to 8 GeV with high resolution. The current status 
of the calorimeter, now in its seventh year of operation, is being presented, as well as details on 
improvements made to the analysis code during the last years. 
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I. The SLAC 5-Factory 

SINCE 1999, the asymmetric B factory at the Stanford Lin- 
ear Accelerator Center (SLAC) is colHding 9.0 GeV elec- 
trons with 3.1 GeV positrons. This total energy of 10.58 GeV 
corresponds to the T(46') resonance. Its decay particles, B^B~ 
and B^B^ pairs, are used to study CP violation and many other 
processes in particle physics. The data are collected by the 
BABAR Detector (Figs. 1 and 2), which is described in detail 
in Ref. [1]. 
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Fig. 2. Schematic drawing of the EAR4R detector. The electromagnetic 
calorimeter is the subsystem shown in purple. 



all pointing close to the interaction point. At the back of 
each crystal two photo diodes and one pre-amplifier card are 
attached (Fig. 3 right). On average, the diodes see about 7,300 
photo-electrons/ MeV. The electronics covers the signal with 
an 18 -bit dynamic range by combining the output of a 10- 
bit ADC with two range bits. This allows the calorimeter to 
measure photon energies from 20 MeV to 8 GeV. The energy 
and position resolution was determined to be the following [2] : 
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Fig. 1. The BABAR detector during the time of construction. 



The crystals are combined into 7x3 modules (except for 
6x3 modules at the backward end and special modules for 
the forward end), then assembled like shown in Figure 4 into 
a barrel and an endcap. Fig. 5 gives a view inside the barrel 
during assembly, while in Fig. 6 the completed barrel is waiting 
for its insertion into the detector. 



II. Electromagnetic Calorimeter Overview 

The electromagnetic calorimeter (EMC) consists of 6580 
CsI(Tl) crystals 16 to 17.5 radiation lengths deep (Fig. 3 left). 
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III. Performance of Hardware 

The operation of the calorimeter is very stable. Out of 
6580 crystals, only one crystal is dead with no hope for any 
recovery. Currently four more crystals are dead, but they might 
be recovered at some time. Fourteen more crystals use only one 
of the two diodes; several more crystals are working incorrectly 
in one energy range, for example at low energy. From time 




Fig. 3. Left: Photograph of a CsI(Tl) crystal lit from the bottom by a light bulb. 
Right: Schematic drawing (not to scale) of a crystal with attached electronics. 




Fig. 4. Cutout drawing of the cylindrical arrangement of the crystals into the 
large barrel and the smaller endcap (left bottom). 



to time an ADC board becomes noisy, which, in the worst 
case, results in the crystals of this board being masked out 
until the board can be replaced during the next access to the 
detector. The electronics is regularly calibrated by determining 
the pedestals and by injecting a known charge into the pre- 
amplifiers to measure the gain and linearity of the system (see 
also Ref. [3]). 

IV. Calibration of Individual Crystals 

The individual response of a crystal to deposited energy, 
namely its overall light yield and its non-uniformity in the light 
output (the dependence on the location of the shower inside 
the crystal), depends on the details at time of manufacture 
and the shape of the crystal. The light output also decreases 
over time due to radiation damage. Each crystal of the BABAR 
calorimeter is therefore regularly calibrated. Two absolute 
energy calibrations are employed for this: The liquid source 
calibration at low energy and the Bhabha calibration at high 



Fig. 5. Photograph of the inside of the barrel calorimeter during construction. 



energy. For intermediate energies the calibration constants are 
interpolations between these two extremes following a function 
linear in the logarithm of the energy. 

A. Liquid Source System 

Whenever a liquid source calibration is performed, a neu- 
tron generator is switched on to emit 14MeV neutrons. The 
generator is surrounded by Fluorinert (FC77), a liquid rich in 
fluorine, and the following chain results in the emission of 
6.13 MeV photons through the decay of ^^N with a half-life 
time of 7 seconds: 
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A system of pipes transports the radioactive liquid past the 
front of the crystals. There the photons enter the crystals and 
are detected with the regular data acquisition system. Figure 7 
shows the spectrum of these photons as seen by a crystal in the 




Fig. 6. Completed barrel calorimeter ready for installation in the detector. 



calorimeter. Calibrations are performed about once a month to 
a statistical uncertainty of <0.5% and a systematic uncertainty 
of about 0.1%. The average loss in light yield over time due 
to radiation damage as measured by the liquid source system 
is documented in Fig. 8. Radiation measurements by RadFETs 
located at the calorimeter indicate that the average radiation 
dose so far is ~0.7kRad for the barrel and ~1.1 kRad for the 
endcap [4] [5]. For more details on the liquid source calibration 
system see Ref. [6]. 




Energy (MeV) 

Fig. 7. Spectrum of 6. 13 MeV photons as detected by a crystal of the calorime- 
ter. The Gaussian functions indicate the contributions from the 6.13 MeV peak 
(far right) and the two escape peaks (middle and left). 
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Fig. 8. Relative loss in light yield over time for the backward barrel (blue), 
forward barrel (red) and endcap (black) as measured by the liquid source 
system. The yellow areas indicate major times without beam. 



B. Bhabha Calibration 

The second absolute energy calibration of individual crystals 
is performed with e+e~ e+e~ events from regular recorded 
data. These calibrations involve crystal energies of 2.5 GeV to 
8GeV, depending on the polar angle due to the asymmetry 
in the beam energies. The calibration requires most crystals to 
have at least 200 direct hits in order to reach a statistical error of 
0.35% for individual crystals. The systematic error is estimated 



to be less than 1% [7]. Calibration constants are currently 
calculated up to once a month, but will soon be obtained more 
frequently once the code is running automatically as part of the 
regular reconstruction system. The constants change over time 
in a way similar to the source calibration constants (Fig. 9). 
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Fig. 9. Evolution of average Bhabha constants over time. 

V. Cluster Calibrations 

Not only calibrations of individual crystals are needed, but 
also calibrations of the clusters, which are groups of adjacent 
crystals in which the full shower energy of a particle is 
deposited. These corrections adjust for shower energy lost at 
the rear of the crystals, gaps between the crystals, and the sides 
of the calorimeter. 

For clusters with energies up to 2 GeV, the calibration is 
obtained from tt^ mesons by correcting the photon energies 
so that the distribution of the invariant mass of two photons 
agrees with the expected tt^ mass distribution. Corrections are 
mostly in the 6% to 8% range. Figure 10 shows a typical 771(77) 
distribution with a clear tt^ peak. An improved version of the 
calibration is currently being tested. 
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Fig. 10. Distribution of the invariant mass of photon pairs after applying the 
cluster calibration. 

For clusters with energies above 2 GeV, the correction factors 
are obtained from single-photon Monte Carlo simulation. Soon 
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VI. Improvements in Reconstruction Software 

A. Position of Cluster Inside Crystals 

When matching a track to a cluster of the calorimeter, the 
position of the cluster in three dimensions has to be known. 
Until recently, this position of the cluster center was always 
located at the front of the crystals. This caused less than 
optimal performance of the matching algorithm, such as when 
a low-energy track, spiraling in the magnetic field of the 
detector, enters the calorimeter at an angle. Now all clusters 
are assumed to be located at a depth of 12.5 cm inside the 
crystals. Figure 1 1 displays the improvement in the matching 
by showing the azimuthal angle difference of the matched 
tracks and clusters before and after this change. Similarly, 
the track-matching efficiency improved, especially at very low 
momentum (Fig. 12). 
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Fig. 1 1 . Distribution of the azimuthal angle difference (in rad) of the cluster 
position and the point where the track intersects with the calorimeter. The plots 
are based on actual data with requirements applied that select a clean set of tt 
mesons. For the top plot the old, for the bottom plot the new cluster position 
algorithm is used. 
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Fig. 12. Track matching efficiency versus momentum for electrons and 
positrons. The plot is based on actual data with requirements that select a clean 
set of electrons and positrons. For the points in red the old, for the points in 
blue the new cluster position algorithm is used. 



data set of the shown distributions is Monte Carlo simulation of 
the decay K*^^. The left plot, obtained without edge 

correction, has a FWHM/2.36 of (45.1 ± 0.7) MeV, while the 
right plot, obtained with edge correction, has a FWHM/2.36 
of (42.0 ± 0.6) MeV, which means that in this case the edge 
correction improved the resolution by 7%. 
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Fig. 13. Ratio of reconstructed over true energy of photons versus crystal ring 
of the barrel. The forward end (where the endcap is attached)_is on the left 
side, the backward end on the right side. The plot is based on BB Monte Carlo 
simulation with measured photon energy ^^cai between 0.5 GeV and 0.8 GeV. 



B. Edge Correction 

If a photon hits the calorimeter at a position close to the 
edge between two crystals, up to 3% of its energy is lost in 
gaps, as can be seen from Figs. 13 and 14. A correction is 
now applied to the energy of each cluster. The effect of this 
so-called "edge correction" on physics analyses can be seen 
from the distributions in Fig. 15. The quantity plotted is AE, 
the difference between the measured energy of a. B meson 
candidate minus the known beam energy. Without measurement 
uncertainties, the peak would be exactly at zero. The underlying 



VII. Additional Studies and Future Goals 

Many decay modes are being used to study the performance 
of the calorimeter, like e+e~ fi/j^j events, radiative Bhabha 
events (e+e~ e+e~7), e+e~ 77, D*^ D^j (photon 
energies from 100 to 400 MeV), Uq (photon energies 

from 50 to 250 MeV). A new cluster calibration is about to 
be implemented, and the Bhabha calibration will soon be 
automated to provide more frequent monitoring and correction 
of the radiation damage at high energies. 
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Fig. 14. Ratio of reconstructed over true energy of photons versus azimuthal 
angle in units of crystals. Due to the modular structure in the calorimeter, the 
crystals are folded over into one single group of three crystals. The dips at 
and 3 are deeper because of the larger gap between the crystal modules. The 
plot is based on BB Monte Carlo simulation with measured photon energy 
^cai between 0.5 GeV and 0.8 GeV. 
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Fig. 15. Distribution of the quantity A^^ for signal Monte Carlo events 
B^ X*+7. The left plot was obtained without, the right plot with the 
edge correction. 



VIII. Conclusion 

The BABAR electromagnetic calorimeter operates very reli- 
ably and delivers very good performance for the experiment. 
The damage to the crystals due to radiation is regularly mea- 
sured and calibrated out. Over time, enhancements were added 
to the reconstruction code, and the tweaking of the calibrations 
continues in order to improve the reconstruction of the detected 
particles and ultimately improve the physics analyses of BABAR. 
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